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Abstract 

Within the concept of the dinuclear system (DNS), a dynamical model is proposed for 
describing the formation of superheavy nuclei in complete fusion reactions by incorporating 
the coupling of the relative motion to the nucleon transfer process. The capture of two heavy 
colliding nuclei, the formation of the compound nucleus and the de-excitation process are cal- 
culated by using an empirical coupled channel model, solving a master equation numerically 
and applying statistical theory, respectively. Evaporation residue excitation functions in cold 
fusion reactions are investigated systematically and compared with available experimental 
data. Maximal production cross sections of superheavy nuclei in cold fusion reactions with 
stable neutron-rich projectiles are obtained. Isotopic trends in the production of the super- 
heavy elements Z=110, 112, 114, 116, 118 and 120 are analyzed systematically. Optimal 
combinations and the corresponding excitation energies are proposed. 

PACS: 25.70.Jj, 24.10.-i, 25.60.Pj 



1 INTRODUCTION 

The synthesis of very heavy (superheavy) nuclei is a very important subject in nuclear physics 
motivated with respect to the island of stability which is predicted theoretically, and has ob- 
tained much experimental progress with fusion-evaporation reactions [H [2]. The existence of the 
superheavy nucleus (SHN) {Z > 106) is due to a strong binding shell effect against the large 
Coulomb repulsion. However, the shell effect will be reduced with increasing excitation energy 
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of the formed compound nucleus. Combinations with a doubly magic nucleus or nearly magic 
nucleus are usually chosen due to the larger reaction Q values. Reactions with ^°*P6 or '^^^Bi 
targets are proposed firstly by Yu. Ts. Oganessian et al. to synthesize SHN [3] . Six new elements 
with Z=107-112 were synthesized in cold fusion reactions for the first time and investigated at GSI 
(Darmstadt, Germany) with the heavy-ion accelerator UNILAC and the separator SHIP [U H]. 
Recently, experiments on the synthesis of element 113 in the ™Zn +^°^ Bi reaction have been 
performed successfully at RIKEN (Tokyo, Japan) [5J. Superheavy elements Z=113-116, 118 were 
synthesized at FLNR in Dubna (Russia) with double magic nucleus ^^Ca bombarding actinide 
nuclei [6j. Reasonable understanding on the formation of SHN in massive fusion reactions is still 
a challenge for theory. 

In accordance with the evolution of two heavy colliding nuclei, the whole process of the com- 
pound nucleus formation and decay is usually divided into three reaction stages, namely the cap- 
ture process of the colliding system to overcome Coulomb barrier, the formation of the compound 
nucleus to pass over the inner fusion barrier as well as the de-excitation of the excited compound 
nucleus against fission. The transmission in the capture process depends on the incident energy 
and relative angular momentum of the colliding nuclei, and is the same as in the fusion of light 
and medium mass systems. The complete fusion of the heavy system after capture in competition 
with quasi-fission is very important in the estimation of the SHN production. At present it is 
still difficult to make an accurate description of the fusion dynamics. After the capture and the 
subsequent evolution to form the compound nucleus, the thermal compound nucleus will decay 
by the emission of light particles and 7-rays against fission. The above three stages will affect the 
formation of evaporation residues observed in laboratories. The evolution of the whole process 
of massive heavy-ion collisions is very complicated at near barrier energies. Most of theoretical 
approaches on the formation of SHN have a similar viewpoint in the description of the capture and 
the de-excitation stages, but there is no consensus on the compound nucleus formation process. 
There are mainly two sorts of models, whether the compound nucleus is formed along the radial 
variable (internuclear distance) or by nucleon transfer at the minimum position of the interaction 
potential after capture of the colliding system. Several transport models have been established to 
understand the fusion mechanism of two heavy colliding nuclei leading to SHN formation, such 
as the macroscopic dynamical model [TJ [8], the fiuctuation-dissipation model [9], the concept of 
nucleon collectivization [TOj and the dinuclear system model [11] . With these models experimental 
data can be reproduced and some new results have been predicted. The models differ from each 
other, and sometimes contradictory physical ideas are used. 
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Further improvements on the mentioned models have to be made. Here we use an improved 
dinuclear system model (DNS), in which the nucleon transfer is coupled with the relative motion 
and the barrier distribution of the colliding system is included. We present a new and extended 
investigation of the production of superheavy nuclei in lead-based cold fusion reactions. For that 
we make use of a formalism describing the nucleon transfer with a set of microscopically derived 
master equations. 

In Sec. 2 we give a description on the DNS model. Calculated results of fusion dynamics and 
SHN production in cold fusion reactions are given in Sec. 3. In Sec. 4 conclusions are discussed. 



2 DINUCLEAR SYSTEM MODEL 

The dinuclear system (DNS) is a molecular configuration of two touching nuclei which keep their 
own individuality [jjj . Such a system has an evolution along two main degrees of freedom: (i) the 
relative motion of the nuclei in the interaction potential to form the DNS and the decay of the 
DNS (quasi- fission process) along the R degree of freedom (internuclear motion), (ii) the transfer 
of nucleons in the mass asymmetry coordinate t] = {Ai — A2)/{Ai + A2) between two nuclei, 
which is a diffusion process of the excited systems leading to the compound nucleus formation. 
Off-diagonal diffusion in the surface {Ai, R) is not considered since we assume the DNS is formed 
at the minimum position of the interaction potential of two colliding nuclei. In this concept, the 
evaporation residue cross section is expressed as a sum over partial waves with angular momentum 
J at the centre-of-mass energy Ecm.i 

(TER{Ec.m) = TT-ET- E (2^^ + J)PcN{Ec.m. J) W^snr (^c.m. , J)- (1) 

Here, T{Ec,m., J) is the transmission probability of the two colliding nuclei overcoming the Coulomb 
potential barrier in the entrance channel to form the DNS. In the same manner as in the nucleon 
collectivization model [10], the transmission probability T is calculated by using the empirical 
coupled channel model, which can reproduce very well available experimental capture cross sec- 
tions [ini [12] ■ PcN is the probability that the system will evolve from a touching configuration 
into the compound nucleus in competition with quasi-fission of the DNS and fission of the heavy 
fragment. The last term is the survival probability of the formed compound nucleus, which can be 
estimated with the statistical evaporation model by considering the competition between neutron 
evaporation and fission [12j. We take the maximal angular momentum as Jmax = 30 since the 
fission barrier of the heavy nucleus disappears at high spin [13]. 



In order to describe the fusion dynamics as a diffusion process in mass asymmetry, tlie analyt- 
ical solution of the Fokker-Planck equation [TT] and the numerical solution of the master equation 
[m [15] have been used, which were also used to treat deep inelastic heavy-ion collisions. Here, 
the fusion probability is obtained by solving a master equation numerically in the potential energy 
surface of the DNS. The time evolution of the distribution function P{Ai, Ei,t) for fragment 1 
with mass number Ai and excitation energy Ei is described by the following master equation 

mm, 

dP{A,, El, t) ^ ^ WA,.A'^{t) [dA,P{A[, E[, t) - dA'PiA,, E„ t) - 

A^^(e(t)) + A^^^(e(t))l p(Ai,i?i,t). (2) 



dt 



Here Wai,a[ is the mean transition probability from the channel {Ai,Ei) to {A[,E[), while dA^ 
denotes the microscopic dimension corresponding to the macroscopic state {Ai,Ei). The sum is 
taken over all possible mass numbers that fragment A[ may take (from to A = Ai + A2), but 
only one nucleon transfer is considered in the model with A[ = Ai±l. The excitation energy Ei is 
the local excitation energy el with respect to fragment Ai, which is determined by the dissipation 
energy from the relative motion and the potential energy of the corresponding DNS and will be 
shown later in Eqs.(8-9). The dissipation energy is described by the parametrization method of 
the classical deflection function [TSl The motion of nucleons in the interacting potential is 
governed by the single-particle Hamiltonian ^121 [T3]: 

H{t) = Ho{t) + V{t) (3) 



with 



K UK 

K,K' qk^Pki 
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Here the indices K,K' {K,K' = 1,2) denote the fragments 1 and 2. The quantities e^j^ and 
UaK,i3K' represent the single particle energies and the interaction matrix elements, respectively. 
The single particle states are defined with respect to the centers of the interacting nuclei and are 
assumed to be orthogonalized in the overlap region. So the annihilation and creation operators 
are dependent on time. The single particle matrix elements are parameterized by 



WaK,/?^/ W = UK,K'{t) I exp 



1 / ^axW 



6. 
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which contains some parameters UK,K'{t) and Ax,K'{t)- The detailed calculation of these param- 
eters and the mean transition probabilities were described in Refs. [HI [T2] . 

The evolution of the DNS along the variable R leads to the quasi-fission of the DNS. The 
quasi-fission rate A'^-^ can be estimated with the one dimensional Kramers formula 



Here the quasi-fission barrier measures the depth of the pocket of the interaction potential. The 



local temperature is given by the Fermi-gas expression O = yE*/a corresponding to the local 
excitation energy e* and level density parameter a = A/12 MeV^^. u^if is the frequency of the 
inverted harmonic oscillator approximating the interaction potential of two nuclei in R around the 
top of the quasi-fission barrier, and u is the frequency of the harmonic oscillator approximating 
the potential in R at the bottom of the pocket. The quantity F denotes the double average width 
of the contributing single-particle states, which determines the friction coefficients: •jn' = ^f^w, 
with yUjj/ being the inertia tensor. Here we use constant values F = 2.8 MeV, huj^if = 2.0 MeV and 
hu = 3.0 MeV for the following reactions. The Kramers formula is derived at the quasi-stationary 
condition of the temperature G(t) < Bqf{Ai, A2). However, the numerical calculation in Ref. [21] 
indicated that Eq.(6) is also available at the condition of G(t) > Bfjf{Ai, A2). In the reactions of 
synthesizing SHN, there is the possibility of the fission of the heavy fragment in the DNS. Since 
the fissility increases with the charge number of the nucleus, the fission of the heavy fragment can 
affect the quasi-fission and fusion when the DNS evolves towards larger mass asymmetry. The 
fission rate A-^*** can also be treated with the one-dimensional Kramers formula EOl 



where ujg,s. and ujf are the frequencies of the oscillators approximating the fission-path potential 
at the ground state and on the top of the fission barrier for nucleus Ai or A2 (larger fragment), 
respectively. Here, we take hug^g, = huf = 1.0 MeV, Fq = 2 MeV. The fission barrier is calculated 
as a sum of a macroscopic part and the shell correction used in Refs. [22]. The fission of the heavy 
fragment is not in favor of the diffusion of the system to light fragment distribution. Therefore, 
it leads to a slightly decrease of the fusion probability (seeing Eq.(17)). 

In the relaxation process of the relative motion, the DNS will be excited due to the dissipation 
of the relative kinetic energy. The excited system opens a valence space Aek in fragment K{K = 
1,2), which has a symmetrical distribution around the Fermi surface. Only the particles in the 



states within the valence space are actively involved in excitation and transfer. The averages on 
these quantities are performed in the valence space: 



^'- = ^J^^'K = e^,9K = j^ (8) 

where e* is the local excitation energy of the DNS, which provides the excitation energy for the 
mean transition probability. There are Nk = Qk^^k valence states and tjik = Nk/'2^ valence 

nucleons in the valence space Ae^, which give the dimension d{mi,m2) = . The 

local excitation energy is defined as 

e* = E,- {U{A,, A2) - U{Ap, At)) . (9) 

Here U{Ai,A2) and U{Ap, At) are the driving potentials of fragments Ai,A2 and fragments 
Ap, At (at the entrance point of the DNS), respectively. The excitation energy E'^ of the composite 
system is converted from the relative kinetic energy loss, which is related to the Coulomb barrier B 
[23] and determined for each initial relative angular momentum J by the parametrization method 
of the classical deflection function [181 HH] • So is coupled with the relative angular momentum. 
The potential energy surface (PES, i.e. the driving potential) of the DNS is given by 

U{Ai, A2, J, R; A, /32, e,, 62) = B{A^) + 5(^2) - [b{A) + Kof (-/)] + V{A^, A2, J, R; A, (32, ^1, 62) 

(10) 

with Ai+yl2 = A. Here B{Ai){i = 1, 2) and B{A) are the negative binding energies of the fragment 
Ai and the compound nucleus A, respectively, in which the shell and the pairing corrections are 
included reasonably. is the rotation energy of the compound nucleus. represent quadrupole 
deformations of the two fragments. 6i denote the angles between the collision orientations and 
the symmetry axes of deformed nuclei. The interaction potential between fragment l{Zi, Ai) and 
2(^2,^2) includes the nuclear. Coulomb and centrifugal parts as 

V{A,, A2, J, R; A, P2, 01, 82) = VW(Ai, A2, R; A, (32, 9,, 92)+Vc{A,, A2, R; A, P2, 9i, ^2)+^^^^^^^, 

2/i rt 

(11) 

where the reduced mass is given by /i = m ■ A1A2/A with the nucleon mass m. The nuclear 
potential is calculated using the double-folding method based on Skyrme interaction force without 
considering the momentum and the spin dependence as 



Vn — Co 



F — F 

Po 



J p?(r)p2(r - K)dr + J pi{r)pl{r - K)dr + F,, J pi(r)p2(r - R)dr\ ,(12) 



with 



F- -f N,-Z,N,-Z, 

J^in,ex — Jin,ex T J in, ex j^^ j^^ i \^'^) 



which is dependent on the nuclear densities and on the orientations of deformed nuclei in the 
colhsion [25j. The parameters Co = SOOMeV ■ fm^, fin = 0.09, fex = -2.59, = 0.42, 
/gj, = 0.54, Po = 0.16 fm~^ are used in the calculation. The Woods-Saxon density distributions 
are expressed for two nuclei as 

^'^"■^ = l + exp[(r-3?i(^^i))/«i]' ^^^^ 

and 

P2(r-R) = — 7- ^ so /aw/ r 

1 + exp[(|r - R| - M2{92))/a2\ 

Here 3?i(6'j) {i = 1,2) are the surface radii of the nuclei with 9ftj(6'j) = + PiY2o{9i)), and 
the spheroidal radii Ri. The parameters represent the surface diffusion coefficients, which are 
taken 0.55 fm in the calculation. The Coulomb potential is obtained by Wong's formula [26] , 
which agrees well with the double-folding procedure. In the actual calculation, the distance R 
between the centers of the two fragments is chosen to be the value which gives the minimum of the 
interaction potential, in which the DNS is considered to be formed. So the PES depends only on 
the mass asymmetry degree of freedom 77, which gives the driving potential of the DNS as shown 
in Fig.l for the reaction ™Zn +^°^ Pb at the tip-tip, the belly-belly and at the fixed (0°, 0°) and 
(90°, 90°) orientations. Here, we should note that the tip-tip orientation is different with (0°,0°). 
We rotate ^ for the fragment with negative quadrupole deformation. However, the orientation 
angle 6i is fixed for all fragments. The same procedure is taken for the belly-belly and (90°, 90°). 
The Businaro-Gallone (E.G.) point marks the maximum position of the driving potential on the 
left side of the initial combination rji. Some averaging over all orientations should be carried 
out in the nucleon transfer process. However, the tip-tip orientation which gives the minimum 
of the PES is in favor of nucleon transfer and is chosen in the calculation. For the reaction 
™Zn Pb, the tip-tip orientation (S/„^=20.98 MeV) has lower inner fusion barrier than the 
belly-belly orientation {Bfus=25.71 MeV). However, the belly-belly orientation appears an obvious 
hump towards symmetric combinations (reducing |?7i|), which is in favor of the compound nucleus 
formation against the quasi-fission. Both of the two factors may affect the values of Pcn (seeing 
Eq.(17)). In Fig. 2 we show the comparison of the formation probability of the compound nucleus 
in the reaction ™Zn-|-^°^P6 as functions of angular momenta (i^c.m. =254.08 MeV, E^j^=12 MeV) 
and incident cm. energies (J=0) at the tip-tip and the belly-belly orientations, respectively. The 
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Figure 1: The driving potential of the DNS for the reaction Pb as a function of the mass 

asymmetry rj at the different orientations. 

effects of the coUision orientations on the fusion cross section were also studied in detail by A. 
Nasirov et al. [27] for deformed combination systems. 

After reaching the time of reaction in the evolution of P{Ai, Ei,t), all those components on the 
left side of the B.G. point as shown in Fig.l (a) contribute to the compound nucleus formation. 
The hindrance in the diffusion process by nucleon transfer to form the compound nucleus is the 
inner fusion barrier Bf^s, which is defined as the difference of the driving potential at the B.G. 
point and at the entrance position. Nucleon transfer to more symmetric fragments will be in favor 
of quasi-fission. The formation probability of the compound nucleus at Coulomb barrier B (here 
the barrier distribution f{B) is considered) and angular momentum J is given by 

Abg 

PcNiE,.m., J,B)= J2 Pi^l^ E,, T^^UEcm., J, B)) . (16) 
Ai=l 

Here the interaction time Tint{Ec,m., J, B) is obtained using the deflection function method [28] . 
We obtain the fusion probability as 

PcNiEc.m., J)= j f{B)PcN{Ec.rr,., J, B)dB , (17) 

where the barrier distribution function is taken in asymmetric Gaussian form [T0| [T2] . So the 
fusion cross section is written as 

oo 

cJfus{E,,^) = 5](2J+ 1)T(E,.^., J)Pcn{E,,^,, J). (18) 

The survival probability of the excited compound nucleus in the cooling process by means of 
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Figure 2: Dependence of the fusion probability on angular momenta and incident cm. energies 
in the reaction '^^Zn +^°^ Pb at the tip-tip and the belly-belly orientations, respectively. 

the neutron evaporation in competition with fission is expressed as following: 

J) = PiE-,,,.,J)U (rjEi^TKEl^).' '''' 

where the E^j^, J are the excitation energy and the spin of the compound nucleus, respectively. 
E* is the excitation energy before evaporating the ith neutron, which has the relation: 

El, = E* - - 2T,, (20) 

with the initial condition E^ = Eqj^. is the separation energy of the ith neutron. The nuclear 
temperature Tj is given by E* = aTf — Ti with the level density parameter a. P{Eqj^, x, J) is the 
realization probability of emitting x neutrons. The widths of neutron evaporation and fission are 
calculated using the statistical model. The details can be found in Ref. [12] . 



3 RESULTS AND DISCUSSION 

3.1 Evaporation residue cross sections 

The evaporation residues observed in laboratories by the consecutive a decay are mainly produced 
by the complete fusion reactions, in which the fusion dynamics and the structure properties of 
the compound nucleus affects their production. Within the framework of the DNS model, we 
calculated the evaporation residue cross sections producing SHN Z=110-113 in cold fusion reactions 
as shown in Fig. 3, and compared them with GSI data for 110-112 [1] and RIKEN results [5] for 

Q 




Figure 3: Comparison of the calculated evaporation residue excitation functions and the experi- 
mental data to synthesize superheavy elements Z=110-113 in cold fusion reactions. 

113. The excitation energy is obtained by E^^ = Ec.m. + Q, where Ec.m. is the incident energy 
in the center-of-mass system. The Q value is given hj Q = AMp + AM^ — AMq, and the 
corresponding mass defects are taken from Ref. for projectile, target and compound nucleus, 
respectively. Usually, neutron-rich projectiles are used to synthesize SHN experimentally, such as 
^"^Ni and '^^Zn, which can enhance the survival probability Wsur in E(l-(1) of the formed compound 
nucleus due to smaller neutron separation energy. The maximal production cross sections from 
Ds to 113 are reduced rapidly because the inner fusion barrier is increasing. Within error bars the 
experimental results can be reproduced very well. There is no other adjustable parameters in the 
calculation. Within the same scheme, we analyzed the evaporation residue excitation functions 
with projectiles ^^Ge, ^"^Se, ^^Kr and ^^Sr to produce superheavy elements Z=114, 116, 118, 120 
in Fig. 4. An upper-limit for the cross section producing 118 was obtained in Berkeley [30] . 

In Fig. 5 we show the comparison of the calculated maximal production cross sections of super- 
heavy elements Z=102-120 in cold fusion reactions by evaporating one neutron with experimental 
data [H H]. The production cross sections decrease rapidly with increasing the charge number 
of the synthesized compound nucleus, such as from 0.2 /i6 for the reaction ^^Ca +^°^ Ph to 1 pb 

in 





E* (MeV) 



Figure 4: The same as in Fig.3, but for projectiles ~^Ge, ^"^Se, ^^Kr and ^^Sr in cold fusion 
reactions to produce superheavy elements Z=114, 116, 118, 120. 
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Figure 5: (Color online) Maximal production cross sections of superheavy elements Z= 102- 120 in 
cold fusion reactions based ^osp^ and "^^^Bi targets with projectile nuclei ^^Ca, ^^Ti, ^^Cr, ^^Fe, 
^'^Ni, ™Zn, ^^Ge, ^^5*6, ^^Kr and ^^Sr, and compared with experimental data. 

for "^^Zn +^°^ Pb, and around even below 0.1 pb for synthesizing Z>113. It seems to be difficult 
to synthesize superheavy elements Z>113 in cold fusion reactions at the present facilities. The 
calculated results are in good agreement with the experimental data. In the DNS concept, the 
inner fusion barrier increases with reducing mass asymmetry, which leads to a decrease of the 
formation probability of the compound nucleus as shown in Fig.6. On the other hand, the quasi- 
fission and the fission of the heavy fragments in the nuclear transfer process become more and 
more important if the mass asymmetry (|?7i|) of the projectile-target combination is decreasing, 
which also reduce the formation probability. There appears a little increase for Z>118, which is 
related to the decreased inner fusion barriers of the three systems. The survival of the thermal 
compound nucleus in the fusion reactions are mainly affected by the neutron evaporation energy, 
the fission barrier and the level density. The survival probability has strong structure effects 
as shown in Fig.6. Accurate calculation of the survival probability is very necessary to obtain 
reasonable evaporation residue cross sections. 
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Figure 6: The fusion and the survival probabihties at J=0 as functions of the charge numbers of 
the compound nuclei with the same combinations as stated in the caption of Fig.5. 

3.2 Isotopic dependence of the production cross sections 

The production of the SHN depends on the isotopic combination of the target and projectile in 
the cold fusion reactions. For example, the maximal cross section is 3.5±f | pb for the reaction 
62 +208 p5 ^269 ^ however 15±^ pb for the reaction ^^Ni Ph ^^71 ^ 
[H EI]. Further investigations on the isotopic trends are very necessary for predicting the optimal 
combinations and the optimal excitation energies (incident energies) to synthesize SHN. In Fig. 7 
we show the calculated isotopic trends in producing superheavy elements Z=110, 112 for the 
reactions "^Ni +^°^ Ph and "^Zn +^°^ Ph (squares with lines), and compare them with the results 
of G.G. Adamian et al. [32] (diamonds and triangles) and the available experimental data [1] 
(circles with error bars). We find that the isotopes 63,64,65 ^^^^ ^"^^^Zn are suitable to synthesize 
superheavy elements 110 and 112, respectively. The isotopes ^'^Ni and ^''Zn have larger production 
cross section, which is consistent with the results of G.G. Adamian et al. But for other isotopes, 
the two methods give slightly different results. For example, our model gives that ™Zn has lager 
cross section to produce elements 110 than the isotope ^^Zn. However, the opposite trend is 
obtained by G.G. Adamian et al. Therefore, it need more accurate description on the three stages 
of the formation of SHN. Further experimental data is also required to examine the theoretical 
models. In the DNS model, the isotopic trends are mainly determined by both the fusion and 
survival probabilities. Of course, the transmission probability of two colliding nuclei can also affect 
the trends since the initial quadrupole deformations depend on the isotopes. When the neutron 
number of the projectile is increasing, the DNS gets more symmetrical and the fusion probability 
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Figure 7: (Color online) Isotopic dependence of the calculated maximal production cross sections 
and the corresponding excitation energies in the synthesis of superheavy elements Z=110, 112 for 
the reactions "^Ni +^°^ Ph and ^Zn +'^^^ Pb, and compared with the results of G.G. Adamian et 
al. [32] and the experimental data [H H]. 

decreases if the DNS does not consist of more stable nuclei due to a higher inner fusion barrier. A 
smaller neutron separation energy and a larger shell correction lead to a larger survival probability. 
The compound nucleus with closed neutron shells has larger shell correction energy and neutron 
separation energy. With the same procedure, we analyzed the dependence of the production cross 
sections on the isotopes Ge and Se to produce the superheavy elements Z=114, 116 shown in Fig.8 
as well as on the isotopes Kr and Sr to synthesize the superheavy elements Z=114, 116 with a 
208 target as shown in Fig.9. It results that the projectiles ''^Ge, "^^Se, ^^Kr and ^7,88^^ ^^^^ 
favorable to synthesize the new superheavy elements Z=114, 116, 118 and 120. The corresponding 
excitation energies are also given in the figures. The compound nuclei with neutron-rich isotopes 
^^Ge, ^'^'^^S'e and 84,86^^., ^^^^ near the sub-closure at N=172. These compound nuclei have larger 
one-neutron separation energies, and the initial combinations smaller mass asymmetries leading 
to smaller evaporation residue cross sections. 
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Figure 8: (Color online) The same as in Fig. 7, but for isotopes Ge and Se to produce superheavy 
elements Z=114, 116. 




Figure 9: (Color online) The same as in Fig. 7, but for isotopes Kr and Sr based ^osp^ target. 



4 CONCLUSIONS 



Within the DNS concept, a dynamical model is worked out for describing the production of 
superheavy residues in the fusion-evaporation reactions, in which the formation of the superheavy 
compound nucleus is described by a master equation which is solved numerically and includes the 
quasi-fission of the DNS and the fission of the heavy fragments in the nucleon transfer process. By 
using the DNS model, the fusion dynamics and the evaporation residue excitation functions in cold 
fusion reactions are investigated systematically. The calculated results are in good agreement with 
available experimental data within error bars. Isotopic trends in the production of superheavy 
elements are analyzed systematically. It is shown that the isotopes 63,64,65 jy^^ ^"^^^Zn, "^^Ge, '^^Se, 
^^Kr and ^'^^^^Sr are favorable to produce the superheavy elements Z=110, 112, 114, 116, 118 and 
120 at the stated excitation energies. 

The physical nature of the synthesis of heavy fissile nuclei in massive fusion reactions is very 
complicated, which involves not only certain quantities which crucially influence the whole process, 
but also the dynamics of the process is important. The coupling of the dynamic deformation and 
the nucleon transfer in the course of overcoming the multi-dimensional PES has to be considered 
in the DNS model. The height of the fission barrier for the heavy or superheavy nuclei should 
be more studied, which is mainly determined by the shell correction energies at the ground state 
and at the saddle point [53] ■ It plays an very important role in the calculation of the survival 
probability. Further work is in progress. 
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